Pressure-Induced Indirect-Direct Bandgap Crossover and 
Bandgap Engineering in Ductile Semiconductor Ag3PO;4 with 
Negative Linear Compressibility 

Yang Lu', Shengcai Zhu'?, Eugene Huang’, Yu He'®, Jiaji Ruan’ & Hao Yan’ 


Silver orthophosphate (Ag3PO4) has aroused widespread interest since the discovery of visible-light- 
driven photocatalysis activity. However, its indirect bandgap nature at ambient becomes an obstacle to 
improve the photocatalysis efficiency and use in other application fields. In this study, the indirect-to- 
direct bandgap crossover was found during the pressure-induced phase transition from cubic to trigonal 
phase in Ag3PO4. The observed bandgap narrowing could be partially retained after pressure was 
released to ambient conditions, accelerating to efficiently utilize the solar spectrum. Moreover, Ag3PO. 
was ductile like metals. Unexpectedly, this ion-conducting and ductile semiconductor photocatalyst 
further exhibited the negative linear compressibility (NLC) at pressure range of 7.5-25 GPa. This work, 
combining experimental and theoretical tools, provides an insight into the modulation mechanism of 
the electronic structure for silver-based semiconductors, and may broaden the applications of Ag3PO14. 
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The direct and indirect bandgap characters of a semiconductor 
play a key role in photo-excited processes, including the separation, 
transfer, and recombination of photo-induced electrons and holes, which 
profoundly affect the physicochemical properties and the performances 
in applications.'? For an indirect bandgap semiconductor, the optical 
absorption coefficients and the recombination are lower than in a direct 
semiconductor because the absorption involves not only the incident 
photons, but also phonons. Interestingly, indirect/direct bandgap 
crossover could balance the absorption and recombination and may 
further observe exotic phenomena, such as valley-based properties.”* 
Recently, the crossover from indirect to direct bandgap was widely 
observed in monolayered materials with quantum confinement effect.” 
8 Applying tensile strain or rising temperature, the crossover can also be 
emerged in multilayered materials due to the decoupling interplay and 
lattice expansion.”"'! Pressure is an effective approach to tune the crystal 
structure, lattice parameters, and electronic structure of materials. !2"!4 
Unfortunately, materials under pressure generally undergo a direct-to- 
indirect and not an indirect-to-direct bandgap transition, attributing to 
compressive effect.'*!¢ 


Silver ion conductors are a versatile class of materials with rich 
mechanical, optical, and transport properties.!’*° Owing to the weak 
interplay of the Ag atoms with the environment, these materials are 
usually ductile under stress and possess excellent ionic conductivity, 
such as the soft semiconductor a-Ag2S and super-ionic conductor a- 
AgI.'868 Recently, the visible-light-driven 
photocatalytic activity of the ion-conducting Ag3PO4 semiconductor for 


the discovery of 


water splitting for green energy has attracted much more interesting in 
the exploiting the silver-based photocatalysts.!??° Under Illumination 
with visible light (wavelength < 500 nm), this novel photocatalyst can 
yield an extremely high quantum efficiency of ~90% towards oxygen 
evolution.'? However, the indirect bandgap nature of the Ag3PO. 
semiconductor at ambient conditions may impede the improvement of 
absorption properties, which in turn lowers the photocatalytic efficiency. 
Furthermore, AgsPO4 should possess good ductility like a-AgoS and 
other exotic properties for new applications. 


In this work, combined with high-pressure experiments and density 
functional theory (DFT) calculations, we found that the indirect-to- 
direct bandgap crossover of Ag3PO4 under pressure couples with a phase 
transition from cubic to trigonal. Experimentally, the bandgap of 
Ag3PO4 narrowed from 2.37 eV to 1.82 eV via compression and the 
shrinking bandgap could be partially retained after releasing the pressure 
to ambient conditions. Besides, this material also exhibits good ductility 
and negative linear compressibility. Via the DFT calculation, we 
confirmed the high-pressure phase atomic structure, the indirect-to- 
direct bandgap crossover, and the bandgap narrowing. Our results will 
shed light on materials with a direct/indirect bandgap engineering under 
pressure and may promotes potential applications in new fields. 


Results 


High purity Ag3PO4 powder (Figure S1) was employed to study the 
ductility, atomic structure evolution, and electronic structure change 
through pressure treatments (Figure la). Upon pressurising, the yellow 
powder could be easily made into a sheet between the diamond culets, 
showing excellent ductility (Figure S2). Then, the pressure effects on the 
optical properties of this sheet-formed material were explored by in situ 
UV-vis absorption spectroscopy. Figure 1b and c show the absorption 
Ag3PO4 
decompression, respectively. As pressure increased from ambient to 50 


spectra of the typical slice under compression and 


GPa, the position of the absorption edge shows red-shifting from 500 nm 


to around 650 nm. Meanwhile, the slope of absorption edge clearly 
decreased when pressure exceeded 8.2 GPa, revealing the degeneration 
of the indirect bandgap and direct bandgap. At a pressure range from 
18.9 GPa to 50 GPa, the slope has a subtly increasing trend (Figure S3), 
which may imply a direct-to-indirect bandgap transition. Upon 
decompression, the changes in position for the absorption edge as well 
as the slope are reversible. 
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Figure 1. Optical properties of ductile AgsPO4 semiconductor. (a) Scheme 
of crystal structure and indirect-to-direct transitions for ductile AgsPO. 
semiconductor by pressure treatment. (b and c) UV-vis absorption spectra of 
AgsPO4 under compression and decompression, respectively. (d) Optical 
bandgap as a function of pressure. The indirect bandgap denoted by red 
(spherical) and black (square) data points, and the direct band gap denoted by 
the blue (star) data points. 


The optical bandgap of Ag3POs as a function of pressure (Figure 1d) 
was derived from the absorption spectra via the Tauc method (Figure 
S4). The region with direct bandgap is marked by gray background 
between 7.5 GPa and 15 GPa. Below the critical pressure of 7.5 GPa, the 
indirect bandgap was 2.37 eV, agreeing with the previous report and 
remained stiff, which corresponds to the cubic phase. "° In the gray region, 
the direct bandgap of Ag3POs4 with Trigonal-I phase jumped up to ~2.5 
eV and then decreased with further increase in pressure. After the direct- 
to-indirect transition, the indirect bandgap of the sample monotonically 
narrowed to 1.82 eV with increasing pressure. Meanwhile, the color of 
the sample changed from yellow to orange during compression, as 
evidenced in the inserted photographs of the Ag3POs slice (Figure 1d). 
Upon decompression, the bandgap hysteretically broadened from 1.82 
eV to 2.26 eV, whilst the color back to yellow. Unexpectedly, the 
bandgap of the quenched sample at ambient condition decreased by 5% 
compared with the starting material, which may lead to an enhancement 
of photocatalytic activity. 


The structural evolution should dominate the variations of 
electronic structure, then the indirect/direct bandgap transitions and 
bandgap tuning under pressure. In this work, the Raman spectroscopy 
synchrotron X-ray diffraction (AD-XRD) 
experiments associated with theoretical calculations were conducted 


and angle-dispersive 
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under pressure to gain the local and global structure information. At 
ambient conditions, Ag3PO4 possesses a body-centered cubic structure 
(P-43n, No. 218) with a lattice parameter of ~6.01 A.2??! In this structure, 
the regular [POs] tetrahedra, constituting the body-centered cubic lattice, 
are isolated by highly distorted [AgO.] tetrahedra and each oxygen atom 
coordinates to one phosphorus atom and three silver atoms. Owing to 
the greater electronegativity of P relative to the Ag atom, the P-O bonds 
are stronger than the covalence-like Ag-O bonds. This then results in a 
prominent Raman peak at frequencies of 908 cm’! for ambient Raman 
spectra (Figure S5a, b), which could be identified as the symmetry 
stretching vibration of [PO4] tetrahedron.***> In high-pressure Raman 
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experiments, the vibration trended to a higher frequency with increasing 
pressure, except the turning range near 7.5 GPa with a red shift of 2 cm” 
' (Figure S5a, b, d, e), which indicated the phase transition of Ag3POs. 
Meanwhile, the prominent Raman peaks gradually broadened after a 
sharp jump around the phase transition point (Figure S5c, f), suggesting 
the slight distortion of the lattice under elevated pressure. During 
decompression to ambient pressure, the broadened Raman peaks could 
not be fully recovered, while the frequency of the [POs] phonon mode 
rapidly reverted to the original position except the hysteresis during the 


high-pressure phase. 
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Figure 2. Pressure-induced structural evolution. (a) Synchrotron X-ray diffraction patterns of AgsPO4 under both compression and decompression (A = 0.6199 


A). (b) Rietveld refinements of the experimental (black fork) XRD data at 1.0 GPa, 9.7 GPa, and 26.5 GPa. The blue lines denote the difference between the 
experiments (black) and refinements (red), and the green vertical lines stand for the refined peak positions. (c, d) Pressure-dependent lattice parameters and molar 
volume of AgsPO,. The relationship between the molar volume and pressure were fitted by the 3'¢-order Birch-Murnaghan (BM3) equation of state (EOS) 
P(V)=3B, f0 +2f;) U+3(B -Df-] where f; =[V,/V)” -1]/2 is the Eulerian strain, B, is the bulk modulus, B, is its pressure derivative, and V and 
V, are the unit cell volume and volume at zero pressure. (e) Variation of compressibility for the Trigonal I-II unit-cell axis as a function of pressure. 


Knowledge of the pressure-induced structural evolution can be gained 
from the synchrotron X-ray diffraction (XRD) spectra. Figure 2a shows 
the XRD patterns of Ag3PO4 under compression up to 50 GPa and 
decompression to ambient pressure. All the Bragg peaks shifted towards 
higher angles when pressure was applied up to 5.2 GPa due to the unit cell 
contraction, suggesting that the cubic structure of the low-pressure phase 
could survive above 5.2 GPa. Three trigonal phases (R3c, No. 161) 
appeared in the pressure range of 9.7-50.7 GPa. Particularly, the cubic and 
Trigonal-I phases coexisted around 7.2 GPa. A peak around the 13.5° 
(marked by the star in Figure 2a) gradually merged with the closed peak 
and disappeared under a pressure of 14.8 GPa, indicating another phase 
(Trigonal-II) transition occurred. Subsequently, the peaks at range of 12°- 
14° continue to merge into one peak when pressure up to 26.5 GPa, 
revealing the Trigonal-III phase emerging. Figure 2b shows the typical 
Rietveld refinements of the XRD patterns at 1.0 GPa, 9.7 GPa, and 26.5 
GPa, respectively (Table S1 and S2). Obviously, the good refinements 
confirmed that Ag3PO4 underwent phase transitions from cubic to trigonal 
structures under high pressure. To our best knowledge, the pressure- 
induced phase transition from cubic to trigonal is, for the first time, 
observed in A3BX4—-type compounds. 


The crystal structure refinement parameters of all diffraction patterns 
are summarized in Figure 2c and d. The lattice parameters and molar 
volume of cubic Ag3PO4 decreased with increasing pressure. Unusually, 
the c-axis of the Trigonal-I and II (or Tri. I and Tri. II) phases exhibits 


= negative linear compressibility with -3.02 ~ -2.87 TPa’! (Figure 2c and e). 


Hence, the Ag3POs, has the potential applications in highly sensitive 


) pressure sensors, optical fibers with high shock resistance, and artificial 
) muscles. The anisotropy of the compressive behavior should be ascribed 
) to the relatively compact stacking along the <001> orientation and the 
) distortion of the [AgOx] polyhedron. After the transition from the 
, Trigonal-II to Trigoanl-III phases around 26 GPa, the lattice parameter of 


c shows a steady, while a continues to decrease at a slower rate. From the 


à volume-pressure relationship, as shown in Figure 2d, the percentage 


decreases of the molar volumes from the cubic to Trigonal-I phase around 
7.5 GPa were calculated to be 3%. In addition, a 3" -order Birch- 
Murnaghan (BM3) equation of state (EOS) was employed to fit the 
experimental data belonging to each structure. Then, the EOS parameters 
were determined to have a bulk modulus of Bo = 101.29 GPa and volume 
of Vo = 217.47 A? for the cubic AgsPOs, Bo = 121.04 GPa and Vo = 210.77 


| A3 for the Trigonal I-II structures, and Bo = 228.96 GPa and Vo = 191.80 
: > A3 for the Trigonal-III structure. For cubic Ag3POs, the bulk modulus is 
_ = very small, suggesting that AgsPOs is very soft like some metals and the 
“© ductile semiconductor a-AgS with a bulk modulus of 59.7 GPa. 

) Therefore, the semiconductor AgsPO4 could be applied into the flexible 


26,36 


devices, such as resistive random access memory. 


After decompression to ambient pressure, the cubic structure of 
Ag3PO4 was recovered with hysteresis viewed from the shifted Bragg 
peaks, which agrees with Raman observation. This means that the 
pressure-induced phase transition of AgPO4 is partially reversible. 
Comparatively, the quenched sample possesses the same lattice parameter 
of 6.01 A with the pristine sample, but with broadened peaks (Figure S6). 
Considering the rapid recovery of the [PO.] tetrahedra from the high- to 
low-pressure phase (see Raman results), the hysteresis and broadening 
should be attributed to the slow movement of the silver atoms to their 
original positions in the cubic unit cell, which is related to the hysteresis 
of bandgap during compression. 
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Figure 3. Mechanism and details of phase transitions. (a) Calculated enthalpy 
of the phase transition at 7.5 GPa. Insets are the pressure-induced structural 
change of the 2x2x2 unit cells viewed from the [111 cubic // [001 ]trigonai direction. (b, 
c, and d) Theoretical bond length of Ag-O, Ag-Ag, and P-O as a function of 
pressure, respectively. 


To better understand the phase transition, theoretical calculations of 
Ag3PO4 were conducted using the Stochastic Surface Walking (SSW) 
method, which has been successful in the TiO2 and FeOOH systems.*7“” 
The pathway sampling was carried out initially using the VASP package 
starting from cubic Ag3PO4 (2 Ag3PO4 units per cell) to generate a 
database of pathways that connect Ag3PO4 with the other phases. In total, 
120 pathways were collected from 500 SSW-crystal steps of pathway 
sampling. While the majority of pathways (80) were found to convert cubic 
Ag3POs4 to Trigonal-I phase, other pathways leading to phases such as 
monoclinic were also identified. The first principles DFT calculations 
based on the variable-cell double-ended surface walking method?! were 
then utilized to verify the energetics of every low-energy pathways, from 
which the lowest energy reaction pathway is obtained. The lowest energy 
pathway leaving the cubic Ag3PO4 phase is the connection between the 
cubic Ag3PO4 and Trigonal-I phase, which is a one-step transition, and the 
other transition pathways leaving the cubic Ag3POs phase are energetically 
much higher. The reaction profile is shown in Figure 3a, starting from a 
cubic Ag3POsg lattice to a Trigonal-I phase lattice. Wefound the phase 
transition is much simpler with only the [PO] polyhedrons rotated. To our 
surprise, the phase transition barrier from the cubic phase to Trigonal-I 
phase is almost zero calculated by the NEB method. This was indeed 
observed in our XRD experiment, where the phase transition takes place 
at 7.5 GPa. Furthermore, the calculated XRD patterns and volume- 
pressure relationship are consistent with the experimental results (Figure 
S7 and S8). 


The phase transition of the Ag3POx from cubic to trigonal phases was 
schematically shown in the insets of Figure 3a, viewed from the [111 ]cubic 
// [001] trigonal Orientation. Apparently, the applied pressure induced the 
change of the [AgO4] polyhedrons, then resulted in the reversed rotation 
with around 10° for alternative [POs] tetrahedrons along the [111 ]cubic // 
[001 ]Trigonai direction. To obtain more details, the calculated bond lengths 
of Ag-O, Ag-Ag, and P-O of Ag3PO4 under pressure are summarized in 
Figure 3b-d. Obviously, all the bond lengths split due to the rearrangement 
of atoms after the cubic-trigonal transition, revealing the symmetry 
breakage of the [AgOx] polyhedron. As pressure increased, the bonds of 


the trigonal phase generally shrunk. Notably, one of the P-O bonds of the 
cubic and trigonal phases narrowed linearly with increasing pressure, 
which is reflected in the Raman vibration (Figure 3d and S5). On the other 
hand, along the [001] trigonal direction, the regular [PO4] tetrahedrons were 
slightly squashed by 1-2% and the [AgOx] polyhedrons underwent a huge 
distortion, resulting in the compact stacking of Ag3POs trigonal structures. 


To gain insight into the substantial correlation of the crystal structure 
with optical properties, the electronic band structures of AgsPO4 with 
different phases were calculated using the DFT method. It is well known 
that the common exchange-correlation functionals, such as the LDA and 
GGA, are severely underestimated in the bandgap. In order to overcome 
this, we used PBEO hybrid exchange-correlation functional to calculate the 
band structure. As shown in Figure 4a, the bandgap at ambient pressure is 
2.42 eV, which is very consistent with 2.37 eV experiment. On another 


hand, the energy band configurations of the cubic and trigonal phases are 
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very different, implying their distinction in optical properties. In detail, the 
indirect bandgap of the cubic phase at O GPa degenerated to a direct 
bandgap at 8 GPa after the phase transition from the cubic to Trigonal-I 
phase. As pressure increased more than 17 GPa, the indirect bandgap 
reappeared in the Trigonal-II and III phases. Figure 4b summarizes the 
pressure-dependent indirect and direct bandgap of AgsPO4. Obviously, the 
energy difference between the direct and indirect bandgap for Ag3PO4 
decreased from ~120 meV to 0 eV when the cubic phase shifted to 
Trigonal-I. Subsequently, the difference gradually increased to ~80 eV at 
44 GPa for the Trigonal-III phase. The indirect/direct bandgap evolution 
trend agrees with the experimental results. 
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Figure 4. Electronic band structure and relation with the coordination of silver. (a) Calculated band structures for AgsPO. with cubic phase (0 GPa) and trigonal 
phases of | (8 GPa), II (17 GPa), and III (44 GPa). (b) Theoretical direct/indirect bandgap of AgsPOs. as a function of pressure. (c) Pressure-induced coordination change 


of silver in [AgO,] polyhedrons. 


The indirect-to-direct bandgap crossover usually occurs in transition 
metal dichalcogenides from bulk to a few-layers especially to monolayer, 
such as MoS, which originated from the quantum confinement effect.*>-7 
Besides, the applied tensile strain and temperature rise can also effectively 
drive the bandgap crossover, from indirect to direct via the common effect 
of crystal lattice expansion.°!! In contrast, pressure generally results in the 
shrinking of the crystal cell and the coupling of neighboring atoms, then a 
direct-to-indirect bandgap transition, such as CdSe and CsPbBr3 under 
pressure.'*!® In our case, the direct-to-indirect bandgap transition of 
trigonal Ag3POs under higher pressure verified the pressure effect on the 
electronic structure. However, it is anomalous that the pressure-induced 
indirect-to-direct bandgap crossover during the phase transition from cubic 
to Trigonal-I phase. 


The yellow Ag3POs can be treated as the black Ag2O added with p- 
block phosphorus element based on the theoretical calculation.!? Hence, 
the size and configuration evolutions of the [AgOx] polyhedron under 


pressure should dominate the electronic band structure and the optical 
properties. Therefore, the prolonged partial Ag-O and Ag-Ag bonds owing 
to the bond length splitting (Figure 3b and c) may primarily degenerate the 
indirect and direct bandgap after the cubic to trigonal phase transition. 
Based on the structural prediction, the coordination number of Ag* and the 
scheme of the [AgOx] group as a function of pressure were obtained and 
are shown in Figure 4c. Apparently, the silver coordination of the 
Trigonal-I phase is equal to that of the cubic phase, but varies the [AgO.] 
structure and breaks its centrosymmetry. When pressure rises up to 44 GPa, 
the silver coordination increases from 4 to 7, accompanying the growing 
complexity of the [AgOx] configuration. Comparatively, the transitions 
between the indirect and direct bandgap for Ag3PO. display a good 
correlation with the Ag* coordination number and [AgOx] configuration. 


In conclusion, the pressure modulations in the crystal structure and 
optical properties of ductile Ag3PO4 semiconductor were systematically 
studied by in situ synchrotron X-ray diffraction and optical absorption 


spectroscopy combined with theoretical calculations. Under compression, 
the indirect bandgap Ag3PO4 semiconductor exhibited a good ductility and 
underwent a reversible phase transitions. The pressure-induced phase 
transition from cubic to trigonal phases in Ag3POs are, for the first time, 
observed in A3BX4-type compounds. The bandgap of Ag3PO4 was 
narrowed from 2.37 eV to 1.82 eV via compression and the shrinking band 
gap could be partially retained after releasing the pressure to ambient 
pressure. In addition, an indirect-to-direct bandgap crossover was found 
during the phase transition from cubic to the first trigonal structures as 
predicted by theoretical calculations. At much higher pressure, the 
transition from direct to indirect bandgap was strongly correlated to the 
weakening in bonds and increased Ag* coordination number. These results 
may provide clues for the understanding of underlying mechanism that 
causes the tuning of the electronic band structure and furnish the 
opportunity in the design of better functional materials. 


Methods 


Diamond anvil cell sample loading. Silver orthophosphate powder 
(Ag3PO4, Alfa Aesar, 99%) was loaded into a Mao-Bell type symmetric 
diamond anvil cell (DAC, 300 um culets) with a ruby chip for pressure 
determinations. T301 stainless steel or Rhenium foil (250 um in thickness) 
was pre-indented to 40 um in thickness and a hole was laser drilled (A = 
1064 nm, Power < 0.3 W) at the indentation center to serve as the gasket 
and sample chamber. Silicone oil was employed as the pressure- 
) transmitting medium (PTM). The pressures were measured with the in situ 
ruby fluorescence technique.*! 


In situ high-pressure structural characterizations. The in situ high- 
pressure angle-dispersive X-ray diffraction (AD-XRD) experiments were 
executed in the Shanghai Synchrotron Radiation Facility (SSRF, BL15U1), 
Advanced Light Source (ALS, BL12.2.2), and Stanford Synchrotron 
Radiation Lightsource (SSRL, BL10-2) with wavelengths of 0.6199, 
0.49594, and 0.6199 A. The resulting two-dimensional patterns were 
integrated using Fit2D software.“ Rietveld refinements were performed 
using the General Structure Analysis System (GSAS) with the user 
interface EXPGUI package.**44 The in situ high-pressure Raman 
spectroscopy experiments were performed on a Raman spectroscopy 
, system in the backscattering configuration with an excitation laser (A = 532 
nm, Renishaw 1000). 


In situ high-pressure optical measurements. The powdered AgsPO4 
sample was pressurised (up to 3 GPa) by DAC without PTM to form a film 
with a thickness of ~10 um. Two pieces of film were peeled off the 
diamond culet and together with a ruby chip were loaded into the DAC 
sample chamber, using silicone oil as PTM. The high-pressure UV-vis 
spectroscopy measurements were carried out by a micro UV-vis 
spectroscopy system (DH-2000-BAL, Ocean Optics). During compression 
and decompression, the Ag3PO4 sample color change was photographed 
by optical microscope (Leica M205C). 


Stochastic Surface Walking (SSW) Global Structure Search. The PES 
of Ag3POs was explored using a DFT based SSW global structure search 
(SSW-DFT). The SSW method*5* is an unbiased, automated approach to 
explore the multidimensional PES of complex systems by taking into 
account the second derivative information“. The efficiency of the method 
for exploring PES has been demonstrated for both aperiodic (molecules 
and clusters*®) and periodic (surfaces**, crystals*!**) systems. The 
algorithm of the SSW global optimization method can be found in 
reference. In detail, the potential energy surface of AgsPOs was carried out 
in a 16-atom (2 Ag3POs units per cell) AgsPOu lattice, using the cubic 
structures of Ag3POs as seeds, which is the most stable phase at ambient 
conditions at 10 GPa. In the SSW search, we generally performed a series 
of parallel runs (4-10 depending on the system) starting from the initial 
guess structures, and up to 300 minima were collected at the first stage, 
from which the most stable configuration was obtained. Next, we verified 
the result with the experiment results. If the most stable configuration of 
the first stage was inconsistent with the experiment, then another 300 


minima were collected. This process was repeated until the configuration 
is consistent with the experiment. Via the ransom search, we found the 
most stable configuration, which agreed well with the experiment (see the 
following section). 


SSW pathway sampling. The SSW algorithm has an automated climbing 
mechanism to drive an atomic configuration from a local energy-minimum 
state to a high-energy configuration along random directions, inheriting 
the idea of the bias-potential driven constrained-Broyden-dimer (BP-CBD) 
method for the transition-state (TS) location." Three steps are involved for 
finding the lowest energy pathway, as described below briefly. (i) 
Pathway collection In SSW pathway sampling, first, we start from one 
single phase (starting phase), and utilize the SSW method to explore all 
the likely phases nearby the phase. A structure selection module is utilized 
to decide whether to accept/refuse once a new minimum is reached. If the 
new phase that is different from the starting phase is identified by the SSW 
crystal method°®, we record/output the IS (i.e. starting phase) and the FS 
(a new phase) of the current SSW step. Then, the program will return to 
the IS by rejecting the new minimum to continue the phase exploration; 
On the other hand, if the new minimum identified by SSW is still the 
starting phase (e.g., the same symmetry but a permutation isomer with 
varied lattice), the program will accept the new isomeric phase and start 
the phase exploration from this phase. We repeated this procedure until a 
certain number of IS/FS pairs were reached. (ii) Pathway screening 
Secondly, we utilized the variable-cell double-ended surface walking 
(DESW) method” to establish the pseudopathway connecting IS to FS for 
all IS/FS pairs®*°°. The approximate barrier is obtained according to 
DESW pseudopathway, where the maximum energy point along the 
pathway is generally a good estimate for the true TS”. By sorting the 
approximate barrier height, we can obtain the candidates for the lowest 
energy pathways. (iii) Lowest energy pathway determination Thirdly, 
the candidate lowest energy pathways are selected to exactly locate the 
true TS by using the DESW TS-search method>’. By sorting the exact 
barrier calculated, the energy difference between the TS and the IS, the 
lowest energy pathways could be finally obtained. All the lowest energy 
pathways will be further confirmed by extrapolating TS towards IS and FS, 
and the TSs need to be validated by phonon spectrum calculation, showing 
one and only one imaginary mode. In this study, we explore the structural 
transformation between the cubic Ag3PO4 phase and the high-pressure 
Ag3POs phases, using the SSW method. Based on the sampled free energy 
landscape, the pathways connecting the initial structure and the final 
structure on the potential energy surface can be determined in an atom-to- 
atom correspondence, and we further export the pathway by the nudged 
elastic band (NEB) method”. 


DFT calculation details. All calculations were performed using the plane- 
wave DFT program, Vienna ab initio simulation package (VASP)*! where 
the electron-ion interaction of Ag, P, and O atoms are represented by the 
projector augmented wave (PAW) scheme and the exchange-correlation 
functional assumes the GGA-PBE type. In the SSW pathway sampling, 
16-atom Ag3POs supercells (both ¢ and Py phase) at 10 GPa were used. In 
order to speed up the PES exploration, we adopted the following parameter 
setting: the kinetic energy plane-wave cutoff 400 eV; and the Monkhorst- 
Pack k-point mesh of (4x4x4) for 16-atom supercells. Once the low- 
energy pathways were located, we switched to high-accuracy calculations 
using a plane-wave cutoff of 800 eV and a denser Monkhorst-Pack k-point 
mesh (8x8x8) for 16-atom supercells. For all the structures, both the lattice 
and atomic positions were fully optimized until the maximal stress 
component was below 0.1 GPa and the maximal force component was 
below 0.01 eV/A. It is well known that GGA severely underestimate 
bandgaps, thus we calculated the density of states using the PBEO hybrid 


exchange-correlation functional. 
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